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Abstract: 


The concept of local dynamical equilibrium In the form of geoatrophlc or 
cycloatrophlc equilibria, which is valid at high altitudes above the solid 
surface, may be adequate to describe the wide range of superrotatlon rates 
observed in planetary atmospheres* It does not address the problem of 
"exactly how this superrotatlon **• is produced and maintained against the 
tendency for friction to oppose differential motions between the atmosphere 
and the underlying planet" (Hide et al,. Nature, I960)* Our analysis, to 
some extent pedagogical, is based on the concept of a thermally driven zonally 
symmetric circulation (e*g. beovy, 1973} Glerasch, 1975)* On planets, with 
rotation axis perpendicular to the orbital plane, the equatorial region is 
heated preferentially* The atmosphere expands, and large scale motions in the 
form of Hadley cells develop which produce a redistribution of mass and 
angular momentum* To provide a basic understanding of the phenomenon simple 
analogies are proposed, where mechanical devices with variable moment of 
inertia are frlctionally coupled to a rotating turntable* Using n three 
dimensional zonally symmetric spectral model and Laplace transformation, the 
time evolution of a fluid leading from corotatlon under globally uniform 
heating to superrotatlon under globally nonuniform heating is simulated* For 
high viscosities atmospheric superrotatlon can be understood In analogy with a 
pirouette* During spin up angular momentum is transferred to the planet* For 
reasonably low viscosities, however, this process is reversed* The Increased 
tendency toward geostrophy, combined with the increase of surface pressure 
toward the poles (due to meridional mass transport), induces the atmosphere to 
eubrotate temporarily at lower altitudes* The resulting viscous shear near 
the surface thus permits angular momentum to flow from the planet into the 
atmosphere where It propagates upwards and, combined with the change in moment 
of Inertia, produces large superrotatlon rates at higher altitudes* 
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INTRODUCTION 


The phenomenon of ntmonphertc supocrotntlon Is ublquttoua In our uolar 
syatom. For Earth, there oxlete evidence suggesting that the rotation rate of 
the thermosphare Is as much ns a factor of P *• 1.4 larger than that of the 
surface <c,g., Xing Hole, 1964). Supecrotation is observed on Mara (e.g. , 
Conrath et al., 1973), Jupiter (e.g., Ingersoll et al., 1979) and Saturn (e.g. 
Alexander, 1962), the most rfrikJ.ng example being that of Venus (e.g., 

Schubert et nl., 1980) where the rototion period of the visible cloud top is 
P « 60 tiroes ahotter than that of the solid body (243 days). Each of these 
atmospheres has a distinctive rigid shell (global average) component of super- 
rotation, irrespective whether the planet rotates In the pcogrode (e.g.. 

Earth) or retrograde (e.g., Venus) direction. 

The difficulty of understanding superrotatlon has recently been 
emphasized in a paper by Hide et al., i960 from which we quote: ”The density 
of the atmosphere decreases rapidly with height and more than 95t of the mass 
of the atmosphere Is contained within the troposphere, which extends from the 
ground to an altitude of ~ 10 km, and the lower stratosphere, Tltese regions 
of the atmosphere rotate faster on average than the underlying solid Earth 
with a typical relative speed of about 10 vn s""^. Exactly how this super- 
rotation (a term used for the generally more rapid relative motions found In 
the very tenons regions at heights around 300 km) Is produced and maintained 
against the tendency for friction to oppose differential motions between the 
atmosphere and the underlying planet is a central problem of dynamical 
meteorology.” 

Numerous theoretical models have been proposed to describe the effect 
(e.g. Schubert and Whitehead, 1969{ Thompson, 1970; Gterasch, 1970; Hide, 

1966, 1971; Gold and Soter, 1971; Gletasch and Stone, 1968; Gleraach, 1975; 
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TiCovy, 1973; the review of RlshbcCht 1972), but none la fully accepted* While 
different procesacs arc probably Important In different regions of the 
atmospheres and In various planetary settings, one feature all share In common 
is differential heating by solar radatlon which has long been known to produce 
pressure variations and zonal wind velocities (e*g. Lorenz, 1967). It Is the 
purpose of this paper to examine the angular momentum budget of this mechanism 
and to provide a mechanistic Interpretation that may contribute to a basic 
understanding of the average or rigid shell component of superrotatlon* The 
Important problem of differential rotation and equatorial superrotatlon will 
be discussed In a subsequent paper (Mayr ct al*, I98U). 

Leovy (1973), Glerasch (1975) and others (e.g. Turlkov and Challkov 
(1971), Kalvay de Rivas (1975)) proposed models for the zonally symmetric 
circulation with superrotatlon resulting from differential heating* It was 
suggested by Leovy and shown by Glerasch that superrotatlon Is maintained in a 
balance between upward transport through the meridional cell and downward 
transport through eddy processes* With a steady state solution of the 
nonlinear energy, continuity and momentum equations, Glerasch described the 
superrotatlon In Che viscous atmosphere of Venus from the planetary surface 
upwards. The zonal wind field thereby rigorously satisfied slip-free and 
stress-free boundary conditions. Due to the neglect of the viscous stress In 
the meridional momentum equation, however, the model Is not entirely self 
consistent in that the meridional winds could not be made to vanish at the 
surface. Our results, which do not suffer from this deficiency nevertheless 
confirm the Importance of the meridional circulation for maintaining 
superrotatlon. 

What remains then to be shown Is how, precisely, the flow Is built up or 
superrotatlon is produced and If avid how the atmosphere receives angular 
momentum In this process? In answering these questions we shall discuss In 
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the first part of this paper simple mechanical analogies that elucidate the 
basic process* In the scconc' part we shall present a spectral model which 
describes in the linear regime the thermally driven zonally symmetric 
atmospheric circulation and superrotatlon* The approach Is taken there to 
study the hypothetical spin up of a planetary atmosphere when differential 
solar heating Is suddenly applied and follow the time evolution until a steady 
state Is reached. Solutions obtained over a wide range of parameters thus 
provide Insight Into the processes that establish and maintain superrotatlon. 
We find that, consistent with the mechanical analogies, two aspects are 
Important: With high viscosity the angular momentum for superrotatlon is 

primarily derived from the change In the moment of Inertia due to mass 
redistribution by the meridional circulation. With reasonable low vlsoclslty, 
however, this mass redistribution and the resulting pressure build-up In the 
polar region at lower altitudes causes the atmosphere to subrotate temporarily 
during the spln-up process so that the atmosphere can pick up the angular 
momentum from the planet giving rise to superrotatlon. 

LOCAL DYNAMICS 

We consider simplified equations to describe the zonally symmetric com- 
ponent of horizontal momentum conservation 
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with U, iionnl velocity 


V, merlcilonnl velocity 

p, naB6 density 

n, kinematic viscosity 

Pi kinetic ptcssure 

aipi rotation rate of the planet 

Qi latitude 

Assuming IVI « lUland with a given pressure flcldi p(r, 0), the solutions of 
(i) and (2) yield in the limit of small but finite viscosity (Mayr et al«, 

198 lb) the equation for gradient flow 

V » -mprcosO + J Wp^^cos^O - |g • (3) 

It reduces to 
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in the geos tropic and cyclostrophic appromlnations respectively. 

Formally, the above solutions represent a local balance between inertial 
forces and pressure gradients. They represent the simplest and perhaps niost 
powerful concept for describing superrotation. On most planets in the solar 
system the rotation axis is perpendicular to the orbital plane, the equatorial 
region is heated preferentially and the pressure tends to increase toward the 
equator < 0) at higher altitudes. The propensity toward superrotation In 
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planotnry atmospheres Is thus explicable In terms of Equations (3) through 
(5). Moreover, the largest superrotation rates are likely to occur at higher 
altitudes away from the viscous constraint by the rotating planet* Tn that 
regime the above solutions are approximately valid and may provide a good 
quantitative description of the phenomenon. An outstanding example is the 
work by Leovy (1973) wViere the large superrotation rate of Venus was 
attributed to cyclostrophlc balance. 

Using the same horizontal pressure gradient, "Iq " “0.3, for different 

planets, we note that the concept of local dynamical equilibrium produces a 

•'’ide range of values ~ ■ 53 (Venus), 0.20 (Earth) and 0.0014 (Jupiter) which 

P 

match the observed average rotation rates within a factor of two* The zonal 
velocity is here expressed in the form of a rotation rate. Am. 

ANGULAR MOMENTUM BUDGET 

The concept of local dynamical equallbrium does not provide a connection 
between superrotation and source mechanism. Moreover, it is restricted to 
that port of the atmosphere where, formally, viscosity is unimportant and thus 
it cannot provide a mechanistic connection between the superrotatlng 
atmosphere and the rotating planet. 

The concept of angular momentum conservation is fundamental for our 
understanding of superrotation. With the primary emphasis on elucidating the 
physics, we shall describe here a scenario of atmospheric superrotation which 
heuristically evolves from fluid dynamical elements that can be readily under- 
stood to a more realistic configuration containing these elements. In support 
of this Interpretation we shall make use of mechanical analogies. 

As reference, a globally uniform atmosphere with mass density Pp(r) Is 
chosen. It is heated with a source, Q(,(r), representing the global average of 
solar input. Due to viscous interaction this atmosphere corotates with the 
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nnguXar velocity of the planet, Wp. W ithout Adding an cxteninl torque to the 
en tire atmosphero-plnnct syatem, w change the heat source to Q(r,0) • 

Q^(r) + AQ(r,0) which In turn changes the density distribution to p(r,G) ■ 
p^^(r) + Ap (r,0). The question we ask Is how AQ maintains an atmospheric 
rotation rate o)(r,0) " “p + Aw(r,0) with a rigid shell component different 
from that o'^ the plonot? 

In t'cality, the rotation rate of the planet also changes during the spin 
up process# For the purpote of discussing superrotation, however, we are 
justified to neglect this oflect since the moment of inertia of the atmosphere 
is much smaller than that of the planet. 

Under conservation of angular momentum, L, the physical conditions for 
uniform and nonuniform heating, representing end states of equilibria, are 
related 

\ (Q) * ^ 


who re 


iL " Lp(Q) - Lp(q„) 


( 7 ) 


represents the total amount of angular momentum that is transferred from the 
atmosphere to the planet (positive) or visa versa (negative) during the spin 
up process# Subscripts a and p refer to atmosphere and planet respectively# 
Assuming Ap « p^, the explicit form of Eq. (6) is 


J " « -Up J' 


Apdij) “ AL# 
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With the definition # “ r cob Odt, the integrations are taken over the entire 
voluroc of the atinosphere (with element dr). Equation (8) states that three 
components contribute to the conservation of the total angular momentum: (a^ 

atmospheric superrotatlon, Aw, (b) changes in the atmospheric moment of 
inertia which are associated with changes In mass density, Ap, and (c) the 
above defined exchange of angular momentum between atmosphere and planet, AL, 
during the transition from one (Qq) to the other (Q) equilibrium state* 

A simple thouglit experiment Illustrates the nature of AL. Starting out 
with an atmosphere character l?,cd by Q^, p^, and w^, corotating with Che 
planet, we lower the heat sonree uniformly over the globe Q(r) ■ QgCr) - 
AQ(r). Under mass conservation the atmosphere contracts, p(r) «• P^Ct) ~ 

Ap(r) and its moment of inertia decreases. During this change, the 
atmosphere tends to jnalntaln its angular momentum and initially there is 
superrotatlon. This motion produces viscous IncernoClon with transfer of 
angular momentum to the planet, until eventually a new static equilibrium Is 
established in which atmosphere and planet corotnte with the angular velocity, 
10 . Of course, for this particular case the solution Is simply 
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to Che planet's moment oi 

: inertia. 

Considering Aw « 0 another 


form of the solution can be obtained from Eq. 8, 
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This moons that the atinospljero transfers an umaunt of angular inomcntuin 
equivalent to the entire change in its moment of Inortlaj thus none remains in 
tlio atmosphere to maintain superrotatlon. 


PtROMKTTK 

V’or a realistic atmosphere with globally nonunlfortn source the conserva- 
tion of angular momentum la complicated, as will bo shown in numerical simula- 
tions of the spin up process. To provide a fundamental understanding of 
superrotation It Is therefore helpful to dlscviss first a simple model which 
proves to l>e valid for very large klncnatlc viscosities. 

The process Is illustrated In Figure 1. f)n a planet where the equatorial 
region ts heated preferentially th;- compressible viscous atmosphere expands. 
Thus large scale motions In the fo.»n of Hadley cells are set up (e.g. Lorenz, 
1967, Leovy, 1973, (Jierasch, 1973) that effectively transport mass from the 
equator toward the poles closer to the rotation axis. As n result the 
atmospheric moment of inertia decreases 


At 
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AQdt < 0 


( 12 ) 


which tends to decrease Che angular niomentura. Under Newton's low such a 
change Is met wldlj inertia and the angular velocity must increase (super- 
rotation) to conserve angular momentum. We assume that during the initial 
phase of the spin up process the atmosphere superrotates all the way down to 
the planetary surface and "reacts” by transferring, through a viscous torque, 
angular momentum to the planet. 

Unlike the case for a uniform source where this moraentura transfer cor- 
responds to the entire, change in moment of Inertia (Eq. 11), only part of it 
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tfl trnnoferred to the planet (AL>0) under nonuniform hontlOR} thus Hupor- 


rotutlon can be sustained, which is restated In mathematical form 


<Aw>I, 


'*<’ ■ -“p / Apd^ -Ah > 0» 


(13) 


Apart from the poleward transport of mass (tenn containtnn; Ap in Eq* 13) 
which is driving superrotation, the meridional winds arc also Important In 
"IsolntlnB*' the superrotntlng atmosphere from the planet. The associated 
Inertial forces (which change sign as one vertically trnnsverscs a Hadley 
cell) are In balance with the viscous forces from the zonal wind field, so 
that in equilibrium auperrotatlon can be maintained with stress free 
boundaries. Solar heating and the resulting meridional circulation thereby 
continually supply the energy that Is expended through viscous dissipation. 
During the spin up process (non-eqnlllbrlum) the inertial forces effectively 
limit the transfer of angular momentum to the planet. 

The imignitude of this momentum transfer is proportional to the viscosity 
and It tins the effect of damping superrotntion. A more viscous atmosphere Is 
tighter coupled to the planet which tends to constrain it Into corotatlon. 

In support of this interpretation a mechanical model la presented which 
Is closely related to the spin up In a pirouette. This analogy Is crude by 
comparison with the dish pan experiments (e.g. Schubert and l^iltehead, 1969) 
providing n fluid dynamic simulation of superrotation. However, it contains 
some of the Important elements to describe the process and has the advantage 
that the physics Is more transparent. 

The discrete phases of the pirouette experiment are illustrated in Figure 
2. Sitting on a rotating turntable, representing the solid planet, is a 
mechanical device referred to as "figure" which represents the atmosphere; Its 
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configuration anil inoniont of inertia (with tnnBsea in ahadoil areas) can be 
phnngod by pulling o string. The angular velocity of tbe figure is schemnt- 
ically shown as a function of Lime. In the first picture frame the string is 
relaxed and no external force is applied. The figure has the maximum moment 
of inertia attoinablo and corotneos with the turntoble, wlvich is the analogy 
of nn atmosphere under globally uniform heating. 

In the second frame the string is pulled. Expending energy against 
gravity, Inertial forces (e.g, centrifugal and Coriolis) and friction, the 
figure utralghtens out, its moment of inertia decreases and it superrotates. 
Locked in this position, the new figure would continue to superrotate indefin- 
itely If there were no frictional Interaction or momentum exchange with the 
turntable (thin solid line In the velocity diagram). Under realistic 
conditions, however, friction causes a slowdown (heavy solid line in the 
velocity il Ingram) which eventually would lend to corotation. To prevent this 
tho ‘ : jt lx» spun up periodically. This is accomplished by relaxing tbe 

string so that the figure assumes its original moment of Inertia. All of the 
potential energy and part of the Initial energy previously expended are 
thereby released. In that position the figure subrotates and picks up angular 
momentum from the turntable. By pulling the string the figure is then spun up 
again, and a new cycle begins. 

The rotation rotes in Figure 2 illustrate the steady state condition 
during nn experiment cycle and indicate that the net result of the procedure 
is superrotation. A simple analysis of the mechanism can readily verify 
this. 

The equation of motion has the form 

(lo,) = - ^ Iio, (U) 
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where I'r) is the time varying moment of inertia^ w Is the difference between 
the angular velocities of the figure (with mass m) and the turntable, and a is 
tlie friction coefficient • Over the experiment cycle, in steady state, the net 
amount of momentum exchange with the turntable must be zero 


o 



^(lu) it 




Itiidt « 0. 


(15) 


and ( 1)2 t’sing the average angular velocities during the two phases of the 
cycle and tp Ij , and t 2 » I 2 being the corresponding time periods 
(t = tj^ + 12 ) and moments of inertia, Equation (15) becomes 


^1 “l 


*^1 + ^2 
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0 . 


(16) 


It foilov’s then Immediately that the average superrotation rate is 
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which is positive, since > u)j ^2 ^ 0 ^nd Ij < The angular vel- 

ocities and ujj 2 sre defined in Figure 2^ The figure superrotates because 
of the cyclic change in the moment of Inertia. Energy is thereby continually 
fed into the system where it is dissipated through friction. The super- 
rotation rate increases inversely proportional to friction. 

Relative to the initial state of corotation (with w^), the subsequent 
state of superrotation has a lower (average) angular momentum. The dlffer- 
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cnce, AL > 0 (Eq. 13) Is trausferrt'd to the turntable during the spin up 
process. It increases with friction and decreases superrototion. 

In summary, we note the analogies with the superrotating atmosphere. The 
oscillatory motions within one cycle of the pirouette experiment simulates the 
Hadley cell. The higher the frequency is in the experiment cycle, the shorter 
is the turnover time in the Hadley cell and the more energy is going through 
the superrotation engine. In the former, mass and momentum exchange take 
place sequentially, while in the latter they occur simultaneously in the upper 
and lower legs of the circulation. The external force on the string repre- 
sents the horizontal pressure gradient due Co differential heating. These 
forces drive Che radial motions of the pirouette experiment and the meridional 
motions in the otmosphere. thereby cone rolling the moments of Inertia in both 
systems. They are primarily balances by Coriolis and centrifugal forces which 
provide the link, with friction and viscous dissipation where energy is 
ultimately dissipated. Thus closing the energy chain we obtain a basic under- 
standing of superrotation which is consistent with the interpretations of 
Leovy, 1973 and Glerasch, 1975. 

One easily sees that the concept of this mechanical analogy can also be 
applied to individual particles of air moving with a Hadley cell. In a fluid, 
at any given time, every phase of this motion is equally represented and thus 
justifies the time averaging procedure as a means of representing the dynamic 
system in steady state. 

DOUBLE PIROUETTE 

The problem with the above analogy is that the entire atmosphere is 
portrayed as a body which can be deformed horiaontally but not vertically. 

Such a picture Is essentially valid when large viscosity provides enough 
rigidity. Under more realistic conditions with reasonably low viscosity. 
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however^ the situation Is more complicated* 

A logical extension of the mechanical model Is the "double pirouette" 
experiment which Is significantly different from the earlier single pirouette 
experiment In that is has an additional degree of freedom for the vertical 
momentum distribution. Instead of one mechnlcal device, two are used on top 
of each other, separated by a friction plate (Figure 3). Pulling the string 
at the upper figure (F^) supplies the energy. As F^ Is stretched and Its 
moment of inertia decreases, the lower figure (Fj_) Is compressed and Its 
moment of inertia Increases. Springs in Fj^ may provide a restoring force due 
to the compressibility of the atmosphere. 

The analytical solution for this thought experiment Is straightforward. 
But depending on the conditions, the results may differ significantly, and a 
detailed discussion would go beyond the scope of this paper. Adopting here a 
simplified description, we consider the corotating (m^) double pirouette with 
moments of Inertia and I^ representing the atmosphere under globally 

uniform heating. The string is pulled during the "energization phase" and the 
moments of inertia change to Ij^; It Is relaxed during the "dissipation 
phase" and the moments of Inertia return back to I^, I^. Subscripts u and 1 
refer to the figures F^ and Fj^, while superscripts e and d Indicate the 
energization and dissipation phases respectively. During the energization 
phase, the time average rotation rates (relative to are w® > 0 and 
< 0 In Fy and Fj^ respectively. 

Conservation of angular momentum with a system of two coupled equations 
corresponding to (15) yields 
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a^, friction coefficient between turntable and pirouette device 
a, friction coefficient between both pirouette devices. 

Since 1° < > 0) and x'r > I'J (w° < 0), it follows from (18) and (19) 

tliat the average rotation rates for both figures are positive 

\ - f <“u > 0 <2'^> 

“l “ T ^“l “l ^ ° 

Exnniing the angular niomentum budget (Gqs. 8 and 13) 

0 < cot 5 w + io, - (0 (I® ~ + I® - 1^) - AL (22) 

o uu 11 o u u 1 1 '' 

produces an interesting result. Suppose we chose I® ■ 1^ and I® » which is 

the case Illustrated in Figure 3. Under this condition the total change in 

the moment of inertia is At 5 (l^ - + I® - I^) «» 0. Knowing from (20) and 

(21) that there is superrotntlon» Eq. 22 then implies that the double 

pirouette device is capable of picking up angu%u- •-..omentumi AL < 0, from the 
turntable during the transition from (uniform) corotation to superrotation. 
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Thl« rofiulc Lb completeXy tllPfercnt from that obtained with the single 
pirmietco experiment where the decrease In moment of Inertia causes 
Buperrotatlon and, additionally, permits angular momentum to be transferred to 
the turntable during the spin up process. Depending on the amount of 

friction, the total change in moment of inertia, AI, and or the (initial) 

exchange of angular momentum with the turntable (planet), AL, may thus con- 
tribute to the angular momentum for superrotation. Under atmospheric 
conditions with reasonably low viscosity it is the latter which is usually 
more important. 

The particular form that was used here to energize the double pirouette 

experiment is not unique. An external force applied to the lower figure would 

load to virtually the same conclusions. 

On the basis of this analysis which is consistent with the results from 
our numerical simulation of atmospheric superrotation in the second part of 
this paper, the following conclusions are drawn: 

Viscosity . In the case of Identically zero viscosity or friction, the 
equations (20) and (21) cannot be applied. But assuming that initially the 
mechanical device (or the atmosphere under uniform heating), somehow, 
uniformly corotates, one can easily show that 5^ > 0 but Wj < 0. The angular 
momenta in regions of super-and sub-rotation add up to zero implying that the 
average change in the rotation rate is also zero. Superrotntion cannot be 
understood without viscosity. To sustain this motion the system must be 
capable of converting kinetic energy back into thermal energy. 

Mass Transport. Without an external torque, conservation of angular 
momentum requires that the moment of inertia must change in order to sustain 
superrotation. This is accomplished through mass redistribution by the 
meridional circulation of the atmosphere induced by solar heating and pressure 
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gradient rt, or by Uw’ osclUntory motions due to extennnl forcing In the 
prlouette exporiraonts. To understand superrocntlon Lt is importtuit thnt the 
change in moment of inortin Is not maintained Btntlcnlly but through transport 
procesBcs. The mass carryin;* meridional circulation of the atmosphere and the 
osclilatory motions of the ini*chanical analogies continually supply the 
energies that are expended by viscous dissipation associated with super- 
rotat Ion. 

Momentum excluinge with the planet . During the transition from corotation 
under globally uniform heating to superrotation under dtfferentinl heating the 
atmosphere exchanges angular momentum with the planet. In simplest form, one 
expects a behavior similar to that of the single pirouette experiment where, 
through viscous Interaction, angular momentum Is transferred to the turntable 
thus slowing down superrotation. However, under more realistic conditions 
with reasonably low viscosity, a secondary circulation develops temporarily at 
lower altitudes and the analogy with the double pirouette experiment is 
appropriate. Initially, during the spin up process of the atmosphere or 
during the energization phase of the mechanical analogy, subrotation prevails 
near the surface, so that the systems can actually receive angular momentum 
from the planet or tVie turntable. 

Analogy with subrotation. The double pirouette experiment represents a 
mechanical model for an atmosphere, preferentially heated at the equator. Our 
simplified description of this analogy, aimed toward a basic understanding of 
superrotation, however, cannot be adopted for subrotation. In reality, the 
changes in the moment of inertia due to the oscillatory motions of the pir~ 
cuette device are continuous and the Coriolis force must be considered. Under 
this condition the analysis Is mors complicated but shows that external com- 
pression against the action of the Coriolis force indeed leads to sub- 
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rotation. This analogy was used to lllustrato the subrotation on the sun 
whose atmosphere Is preEerentlalXy heated at the poles (Schatten et al«, 

1980). Analogous to the fluid <e.g. Eq. 3), the mechanical experiment is not 
symmetrical for large values of super- and sub-rotation when the centrifugal 
force IxicomCB important. 

FLUID PyNAMIC MODEL 

In the framework of a spectral model, we shall concentrate here on the 
lowest order harmonic or rigid shell component of supercotation which Is 
unique in that it alone pocesses a global angular momentum; the angular 
momenta of the higher order vector spherical harmonics are Identically zero. 

A subsequent paper will discuss the problem of differential rotation and the 
acceleration of the equatorial region (Hide, 1970) by extending the analysis 
to higher order modes under consideration of mode coupling and viscous flow 
(Mayr et al., 1981a). 

In the following sections the model is presented, and boundary conditions 
employed at the planetary surface are discussed. Steady state solutions are 
then considered followed by a presentation and discussion of time evolution 
simulations. The thrust of this study is toward a general understanding of 
superrotation rather than to provide quantitative results for a specific 
planetary atmosphere. 

To describe atmospheric superrotation, a three dimensional spectral model 
(Mayr and Volland, 1972; Harris and Mayr, 1975; Mayr et al., 1978, referred to 
as MH) is adopted whose properties are summarized here: (1) Energy advectlon 

and convection associated with large scale circulation are considered. (2) 
Heat conduction and dissipative processes due to ion drag and viscosity In 
both the zonal and meridional force balances are Included. (3) Vertical and 
horizontal mass transport including momentum transfer between the major 
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spaclos o£ the atmosphere are considered In self consistent form* (A) 
Perturbation theory la applied relative to n globally uniform atmosphere which 
corotaterj with the underlying planet* The physical quantities including 
temperature » densities and winds are represented In an expansion of spherical 
harmonics and time functions such as Fourier hurmonlcs* which effectively 
separate latitude and time dependences from altitude variations* (5) The 
equations of energy » mass and momentum conservation are integrated from the 
planetary surface up to the exosphere. At the lower boundary, ttie vertical, 
meridional and zonal velocities relative to the rotating planet are set to 
zero, ond adiabatic conditions arc assumed to prevail* For the upper boundary 
conditions, the properties of molecular heat conduction and viscosity require 
that the vertical gradients of temperature and horizontal winds vanish* 

In the tonally symmetric circulation (m«0) we consider the ’’fundamental 
mode" of superrotation 


^^2 

AT, 


X, (r) P? (0) - X, (r) 1/2(3 sln^O - 1) 




V 2 - U2<r) (0) - jg U2(r) sin Q cos 0 

tij - Vj(r) (0) ^ AurcosO 


(23) 


(2A) 

(25) 


^ ^ 

where AQ2 is the external heat source, and AT^, Ap^, W2, V2» Uj are the 

resulting perturbations in temperature, density and velocity components in the 

vertical, meridional and zonal directions respectively* Standard notations 

0 ^0 +0 

are used for the scalar- (P2) vector- (82* Cj) spherical harmonics. Due 
to the orthogonality of spherical harmonics, the contributions to the global 
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angular roomenturo budget (Eq« S) from the zonal wind field and the maia 


redistribution (or changes In the moment of Inertia) arc 

/“o ^ ““ ° ''n'S “l 

"/% ''l «=“ 


(26) 


AoircosOC^dT 



r^cos^ 0 dT £ p° | (2pg- P°) dx 


(27) 


Apart from the corotating globally uniform atmosphere » only the fundamental 
mode or the rigid shell component of rotation (angular velocity Independent of 
latitude) has non~vanlshing angular momenta due to Am and ^P 2 * higher 
order modes which represent differential rotation (variable angular velocity 
as a function of latitude) the global angular momenta are Identically zero at 
each height. 

Wc emphasize again that In this paper no attempt Is made to simulate a 
particular scenario of atmospheric (rigid shell) superrotation which would 
require a thorough analysis of the energy sources (due to e.g« solar 
radiation, auroral particles and Joule heating) and external forces (due to 
e*g. electric fields of solar wlnd-magnetospherlc origin) chat are 
important. For Illustrative purposes we chose a global average atmosphere 
with Tq, pg which corotates with the planetary rotation rate Wp representing 
the conditions on Venus* To drive superrotation, a specific heat 
source AQ 2 Is adopted which peaks at an altitude of 50 km. With these Inputs 
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fixed • the height independent eddy diffusion coeffecient K and the 
corresponding coefficient of eddy viscosity arc varied to examine the effects 
on atmospheric supcrrotatlon. 

BOUNDARY CONDITIONS 

Disregarding the complexities due to corop08ltlon» eight boundary 
conditions are required to solve the ordinary differential equations that 
describe the vertical structure of temperature, density and velocity fields* 

Of these, seven conditions are firmly established: (a) the vanishing vertical 

and horizontal (zonal and meridional) velocities at the lower boundary, (b) 
the requirement that the energy equation is satisfied under neglect of heat 
conduction at the lower boundary and (c) the assumption that the heat 
conduct ion~and momentum-fluxes vanish at the upper boundary* 

In steady state, one must further require for the rigid shell component 
of superrotation that the viscous torque or the momentum transfer vanishes at 
the planetary surface r^ 


(|~)p “ 0; « 0) (28) 

o 

In Implementing the stress free boundary condition (28) it is important to 
realize that mass redistribution (and the resulting change in moment of 
intertla) Is one of the three elements that balance the steady state angular 
momentum budget (see Eq* 8)* This mass redistribution also produces a 
pressure variation at the surface, Ap^, which, a priori, is not known* Also 
contributing to the momentum budget is the atmospheric momentum exchange with 
the planet, AL, during the spin, up process which is related to the surface 
stress. We require, therefore, chat the surface pressure or density 
amplitudes associated with are chosen such Chat the stress free boundary 
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""Sc. £s s*t£sfi*i in st«*iT st*t# 


■•*': r 




?r r. 




ns 


Ir; -^ic c3B?uc*r node* *n iteritire prc^:edure Is us«d vitH M.*%*ry 

se«r;*', z'j enS'izze c^sadltloo C2?^. 

Is suns out that for a slovlr rjsacing plaivet such as Venus this 


riscorous approach to the boundary value prablea produces results reasonably 
olose to those obtained with the (redundant) assucption of hydrostatic 
equillbrlua. This indicates that without constraint the syscea tends to find 
its natural state and in turn justifies using the hydrostatic condition for 
the tloe dependent solution. 

The concept of Influencing the angular ooaentUB budget by varying the 
surface pressure distribution proves valuable in still another way. Thus a 
computer experiment can be performed in which the surface pressure 
distribution Is chosen so that the atmospheric angular momentum is invariant 
during the transition from corotation to superrotatton 


(ip ) + AL * 0. (30) 

^ o 

STEADY STATE 

For a fixed heat source* A which peaks at the equator* Figure 4 shows 
computer results of rigid shell superrotation in steady state ( ■ 0)* 

different values of kinematic viscosity corresponding to eddy diffusion 
coefficients of K » 10*^ and 10^ are chosen. Solid lines represent realistic 
conditions for which the viscous stress or the vertical gradient of the zonal 
velocity Is forced to zero at the planetary surface (Eq. 29). Under this 
condition* rigid shell superrotation prevails for both vlscositleH throughout 
the atmosphere. 
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In analyzing rha angular momentum budget we evaluate the Integrala of 
siiperrotatlon and tnaea redistribution from Eq* 8, or Eqa* 26 and 27 


S Aoj(r)p^r^co8^0dT 

(30 

P2(0)r^co8^0dT 

(32) 


where Ah • M-S. The ratios 8 ^ are given In Figure 4 and reveal that In the 

case of large viscosity the atmosphere has lost angular nonentum to the 
planet AL > Oj while for the smaller value of viscosity the atmosphere has 
gained angular momentum from the planet Ah < 0 during the spin up process* 
Retaining the Input parameters, we show for comparison a second set of 
solutions (dashed lines) which are artificial and diagnostic In nature* Here 
the surface pressure amplitude Is adjusted (Eq* 30) so that, under 
differential heating, Ah»0» For large viscosity superrotation prevails all the 
way down to the surface. In that state the velocities are consistently larger 

than those under the stress free boundary condition (solid line), and a 

d Alt) 

surface torque ^ ^ tends to slow down the excess superrotation while 

transferring angular momentum to the planet* Obviously, such a condition 

cannot prevail In steady state. Initially, during the transition from 

corotation under uniform heating to superrotation under nonuniform heating the 

atmosphere tends to conserve Its angular momentum* Solutions for AL"0 are 

therefore Interpreted as heuristic representations of the transient state 

which, combined with the true steady state solutions “ 0» may serve to 

o 

Illustrate the mechanics of the spin up process. 

For small viscosity the situation is reversed* With AL ■ 0 the 
velocities are now consistently lower than those under the stress free 



boundary condition* Tbo ntmospbere actually subcotatoa (minus) at lower 

altitudes. In this state a negative surface torque (1^)^ < 0 Induces the 

o 

atmosphere to pick up angular momentum from the planet which contributes to 
the larger superrotation* 


TIME EVOLUTION 

Steady state solutions contain no information on causality and therefore 
cannot afford a complete understanding* An example is the Hadl-$y cell: Mass 

transport In the tipper leg of the circulation away from the source and in the 
lower leg of the circulation toward the source are in balance* But in the 
transient process of establishing the circulation^ it is precisely the 
Inbalance of mass flow which is building up the pressure in the polar region 
at lower altitudes, and It Is that pressure force which keeps the atmosphere 
In motion* 

We shall find that in principle the situation is similar for 
superrotation. Tttere is, however, an important additional complication in 
that the atmosphere is coupled to a planet with virtually unlimited angular 
momentum. Wlille the atmosphere's mass budget that governs the Hadley cell 
formation is closed, Its angular momentum budget which governs superrotation 
Is not. 

In describing the spin up process we assume that up to the time t ■ 0 
the atmosphere is heated with a globally uniform source, ” 0, and 
corctates with the planet. Am ■ 0* A globally nonuniform source is introduced 
as a step function 


AQ^ (t) - 0} t < 0 
A (t) ■ const*; t > 0. 


(33) 
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Using Lnplacc tcnnsforroations, the time dopondences of Che heiic source (33) 
and Cho physical variables X(c) (Equ. 23, 24, 25) are described in the form 

00 

4 q2 Cf) “ / A Q 2 " T 

o 

and 

(f) x(f) ■ ^ X(t) dt, (35) 

respectively* With the time dependence f being a variable time constant, 

the NH model Is used to compute the transfer function x(f) for unit heat 
input. The integral equation (35) Is solved with the Gauss-Legendre technique 
of Bellman ct al, 1966. Using six Gausslon points, the calculations are 
carried out over various timescales to provide the necessary resolution* 

For the same input values that produced the results of Figure 4, the time 
evolution of atmospheric superrotation is simulated and is given In Figures 5, 
6, and 7. For large viscosity the atmosphere superrotates at all altitudes 
throughout the spin up process* Immediately above the surface (2 km altitude) 
a transient velocity maximum developcs. In time, this velocity returns to a 
small value near zero. The corresponding (transient) surface stress has the 
effect of decelerating and accelerating the atmosphere and the underlying 
planet respectively* That momentum transfer is also reflected in the steady 
state superrotation rates throughout the atmosphere which are lower than those 
during the initial phase of the spin up process. This confirms our earlier 
interpretation of the solutions in Figure 4* Judging from the time evolution 
of the height integrated superrotation rate we find that the atmosphere spins 
up like a solid body and transfers angular momentum to the planet, in analogy 
With the single pirouette experiment (Figure 2). 

For the smaller viscosity the results arc entirely different (Figure 
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roCiit tonal nnergy comes from the planet ^ a scale analysis yields 

V T K A y. «F I 

S ■ T •* a W 

3 Atk) 

where the symbols K» anti t stand for the viscosity, the vertical velocity 
gradient near the surface during the spin up process and the spin up time 
respectively. A decrease In viscosity is overcompensated by an Increase In 
the. gradient and spin up time so that the atmosphere can receive an 
Increasingly larger amount of angular momentum* The meridional circulation 
transports angular momentum upwards and Its time response to the solar 
radiative Input Is almost Instantaneous. By comparison, the viscous flow of 
angular momentum from the planet Is slow and slows further down as the 
viscosity decreases. Thus with decreasing viscosity, larger velocity shears 
must be built up and a longer spin up time Is required before the viscous flow 
can match the convective momentum transport. 

The dependence of atmospheric superrotation on viscosity Is Illustrated 
in more detail with Figure 8. For a fixed heat source and a wide range of 
eddy diffusion (viscosity) coefficients, the velocities are shown at various 
attitudes. To characterize the angular moraen»:u>Q budget, the ratio between 
mass and rotational momenta, R, is presented. We identify two regimes: (a) 

For lower viscosities and R < 1 (AL < 0) the atmosphere receives an 
increasingly larger amount of angular momentum as K decreases. The velocity 
Is Inversely proportional to viscosity, (b) At the higher viscosities K > 2 x 
iO® the atmosphere is losing angular momentum to the planet and the velocity 
decrease Is more precipitous. For AL > 0 and AL < 0 the analogies with the 
single and double pirouette experiments are appropriate, respectively. 

COWCLUStON 

The spin up process lending to superrotation Is discussed for the linear 
regime of the thermally driven tonally symmetric circulation. Depending on 
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the state of turbulence the atmosphere can go through two different modes. 

For large viscosity (K > 2x10®, In this study) It behaves almost like a solid 
body, tightly coupled to the planet. The superrotation rates are small and 
the spin up times are short. The angular momentum of superrotation Is 
"driven" by the decrease In the moment of Inertia due to mass redistribution 
by the meridional circulation. During the spin up time the atmosphere 
transfers angular momentum to the planet. This process has an analogy In the 
single pirouette experiment. 

For lower values of viscosity there Is an important additional degree of 
freedom In the vertical momentum distribution, and the analogy with the double 
pirouette experiment Is appropriate. Mass can more freely return from the 
poles In the lower leg of the meridional circulation. The change In the 
moment of inertia Is therefore not as large. This return flow combines with 
the increased tendency toward geostrophy to cause subrotation In the lower 
region of the atmosphere during the initial phase of the spin up process. 
Through viscous shear the planet applies a torque to the atmosphere which 
propagates up and in time supplies the angular momentum for large super- 
rotation consistent with a suggestion by Leovy (1973). Tlie spin up times in 
this process are long. Below a critical altitude level which Is about 40km 
for K“10^ the momentum flow during spin up is apparently carried by upward 
propagating gravity and viscosity waves. For a critical viscosity of about 2 
X 10®, in this analysis, the atmosphere does not exchange angular momentum 
with the planet. 

It is again emphasized that no attempt was made here to describe 
differential rotation and equatorial Jets, such as observed on Jupiter and 
Saturn. By expanding the spectral model to include higher order modes, the 
concept of mode coupling Is used to describe differential rotation In an 
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atmosphere with zoaal and meridional viscous stresses (Mayr et al*^ 198ia). 
Momentum is thereby allowed to cascade from lower order modes (which "absorb" 
most of the solar energy) to higher order modes » and the feedback on the rigid 
shell component of superrotation Is considered. At altitudes above the center 
of the meridional circulation, where the specific heat Input is relatively 
large, angular momentum Is effectively advected poleward, except near the 
equator where the meridional winds and the Coriolis force vanish and where the 
horizontal momentum transport must be carried by horizontal diffusion (viscous 
flow). This causes the atmospheric rotation rate to decrease at low latitudes 
and Increase at high latitudes, relative to the rigid shell component of 
superrotatLon; Inducing a tencency for the equatorial region to subrotate. 
However, below the center of the meridional circulation and below the maximum 
energy Input, energy and angular momentum are advected from high to low 
latitudes. Thus the equatorial region can be accelerated to the extent that 
positive jets develop. Again, near the equator, angular momentum imist be 
carried by horizontal diffusion; the equatorial minima In the angular 
velocities of Jupiter and Saturn being evidence fur down gradient momentum 
flow toward the equator* 

In this model the tropospheres of Jupiter and Saturn are assumed to be 
convecclvely unstable with a superadlabatlc temperature lapse rate accounting 
for global scale energy diffusion toward higher altitudes where radiation to 
space becomes important. Under this condition, the upward motions in the 
direct and Indirect atmospheric circulations Imposed by Insolation can then 
supply dynamic energy for "localized" equatorial heating which In turn drives 
equatorial superrotation. Preferential solar heating at low latitudes Induces 
the ordered structure In the atmospheric motions. The efficiency of this 
process Is determined by the stability of the troposphere, the energy 
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conversion In the Interior and the radiative loss from the upper troposphere 
and Btratoaphere. Differences In these global properties may be responsible 
for the large differences between the zonal velocity fields on Jupiter and 
Saturn. Assuming zonal symmetry and thermal forcing by solar differential 
heating, we believe that three factors are important In understanding the 
atmosphere dynamics of these planets; (1) Their large size and rotation rate 
produce through mode coupling a multiple Ferrel-Thomson type circulation with 
upward motions at low latitudes. (2) Energy from the Interior Is diffusing 
radially outwards and maintains a superadlabatlc lapse rate which provides the 
condition for trapping some energy that Is advected toward the equator In the 
lower leg of the meridional circulation. (3) Tlie eddy diffuslvity Is suf- 
ficiently low to require large temperature and velocity gradients for the 
energy and momentum balances, the atmosphere of Saturn being perhaps less 
turbulent than that of Jupiter. 
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Figure Captions 

Figure 1: Planet with rotation axis penpendlcular to the orbital plane around 


the Sun* Schematic diagram Illustrating the energy deposition In 
the equatorial regloni the Hadley cell and the resulting mass relo- 
cation toward the poles. Under the tendency to conserve angular 
momentum the dynamically induced decrease In the moment of inertia 
produces superrotation. 

Figure 2: Single pirouette experiment as analogy for an atmosphere with high 

viscosity. A mechanical device (figure), Its moment if Inertia 
being cyclically changed by an external force (representing solar 
differential heating apd pressure gradient), Is frlctlonally 
coupled to a rotating turntable. Averaged over an ^'experiment 
cycle" the figure superrotates receiving the angular momentum from 
the (average; decrease in moment of Inertia. During the spin up 
process, angular momentum, AL, is frlctlonally transferred to the 
turntable which slows down superrotation. 

Figure 3: Double pirouette experiment as an analogy for an atmosphere with 

(reasonably) low viscosity. Compared to Figure 2 a second figure 
provides the added degree of freedom for the vertical momentum 
distribution. By alternatlngly changing the moments of Inertia of 
both figures they superrotate In the time average* But this can 
now be accomplished without changing the total moment of Inertia, 
thus the device Is capable of picking up the angular momentum for 
superrotation from the turntable. 

Figure A: Computed vertical distribution for the rigid shell component of 

"superrotation", UcosQ. Velocities are zero at the surface. In 
solid lines steady state solutions are presented with stress free 
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boundary conditions* In dashed lines an artificial condition is 
shovm for which the angular momentum of the atmosphere is 
conserved, AL ■ 0* It is interpreted as a heuristic representation 
of the transient state during spin up* The results are given for 
two viscosities corresponding to eddy diffusion coefficients of K * 
10® (Figure 4a) and K - 10^ (Figure 4b). plus and minus signs 
Indicate super- and sub-rotation respectively* 

Figure 5*. Simulation of the time evolution leading to superrotation for K • 
10®. The latitude dependent component of solar heating which 
causes superrotation is turned on abruptly as Indicated in the 
dashed line* Note that the atmosphere super-rotates throughout the 
spin up, with the positive velocity above the surface providing the 
viscous stress for momentum transfer to the planet during that 
process* 

Figure 6: Results for the same heat source as in Figure 5 except that K <■ 

10^* The atmosphere at lower altitudes subrotates during the 
initial phase of the spin up* The negative velocity above the 
surface is picking up a torque form the faster rotation planet 
which provides most of the momentum for superrotation* Note that 
in both results (Figures 5 and 6) superrotatlon prevails all the 
way down to the planetary surface in steady state. With low 
viscosity, the spin up times are longer and the superrotation rates 
are larger compared to the corresponding values for large 
viscosity. 

Figure 7: The results from Figure 3 are shown for mass flow, P^U, at various 

(approximate) times during the spin up process* Note that 
initially, about 2 months after energization, the angular momenta 
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In regions of super- and sub-rotation alsK>st balance each other, 
thus justifying our interpretation of the AL * 0 solutions iu 
Figure 4. 

Figure 8: Dependence of the steady state superrotation rate on viscosity* 

Also shown is the ratio R between the total angular aoBents due to 
changes in the atmospheric monent of Inertia end rotation rate* 

For R » 1 supercotatlon Is primarily driven by the mss 
redistribution (Figure 4a and S) with angular noaentum being 
transferred to the planet during spin up. For R « 1 superrotation 
is primarily driven by the angular momentum picked up from the 
planet during spin up (Figures 4b, 6 and 7). For R > 1 (AL > 0) an 
analogy with the single pirouette can be Mde, while 
for R ^ 1 (AL ^ 0) the analogy with the ante complicated double 
pirouette must be invoked* For the critical value 

a 

of K * 2 X 10 the atmosphere does not exchange angular momentum 
with the planet. 
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